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ABSTRACT 


This thesis examines spectrometric oil analysis data from the A/E35U-3 and the 
Portable Wear Metal Analyzer (PWMA). Data from the two instruments is compared 
using simple linear regression analysis to obtain equations that can be used to convert 
the A,E35U-3 values to PWMA values. The current JOAP Laboratory Manual 1s 
based on A/E35U-3 measurement values and is used as an aid in determining the 
condition of lubricated engine components. The lack of a counterpart manual for the 
newly introduced PWMA makes it mandatory that values based on the A/E35U-3 have 
a conversion to values based on the PWMA. Through transformation equations to be 
developed in this paper, it will be possible for the PWMA user to utilize the current 


laboratory manual until a manual based on PWMA values 1s released. 
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I. INTRODUCTION 


Preventing failures in high performance aircraft engines and reciprocating engines 
due to foreseeable problems is a common sense part of the armed services mandatory 
safety programs. These programs call for regular maintenance to be done to decrease 
the opportunities for equipment failure. Through various programs it was found that 
some failure types can be anticipated by using spectrometers to trace wearmetal 
contamination levels over time in used oil samples. The actions taken based on the 
analvsis of data resulted in savings from reduced failures and decreased maintenance on 
items that were not in need of scheduled periodic maintenance. 

The main spectrometric device currently in use by the armed services is the 
A/E35U-3. It 1s a bulky, heavy instrument that is not easily transportable to 
deployment locations. A newer device, the Portable Wear Metal Analyzer, has been 
designed to fill the need of a ‘deployable’ oil analyzer. Its production model 
distribution is scheduled to commence in the fall of 1987. The nature of the two 
instruments causes them to produce differing measurements of the contaminant 
wearmetals in oil samples. The differences are primarily due to their differing physics 
and engineering design. 

The measurement differences are important to those units deployed with a 
Portable Wear Metal Analyzer. Users of either instrument rely on the Joint Oil 
Analysis Program Laboratory Manual [Ref. 1] which is based on A/E35U-3 
measurements. The oil analyst uses the manual for three primary purposes: to 
compare Wearmetal values obtained with the wearmetal evaluation criteria and trend 
tables; to obtain information from the table providing decision making guidance; and 
to obtain supplemental information providing additional wearmetal diagnostic guidance 
for each type of equipment evaluated. Therefore to use the manual properly, the 
analyst must have the same type values as the technical manual. Through 
transformation equations to be produced in this paper, the Portable Wear Metal 
Analyzer user will be able to convert technical manual values to those that are 
consistent with it. As usage of the new instrument increases, the benefits of including 
separate wearmetal evaluation criteria and trend tables based on the Portable Wear 


Metal Analyzer measurements will become increasingly apparent. 


An analytical development of transformation equations between the Portable 
Wear Metal Analyzer and the A/E35U-3 will be presented in this paper. It will be 
supported by a background of the oil analysis program, descriptions of the two oil 
analysis instruments, and a discussion of the data available. The methodology of the 
analysis and a description of the models used to arrive at the final conclusions of what 
transformation equations are required to translate technical manuals from A/E35U-3 


values to Portable Wear Metal Analyzer values will be given. 


Ii. OIL ANALYSIS PROGRAM 


A. BACKGROUND 

For several years, the Spectrometric Oil Analysis Program (SOAP) has been in 
use by all three military services. The Navy started a trial program in 1956 to 
determine if spectrometric analysis of oil samples could be used to predict engine 
failures. The success of the trial program led to the Navy Oil Analysis Program 
(NOAP). The inter-service importance of oil analysis was seen and culminated in the 
1976 merger of the Army, Navy, and Air Force programs into the Joint Oil Analysis 
Program (JOAP). 

Of many statements made on the purposes of the Joint Oul Analysis Program, 
Thomas Menard provides one that is simple and concise. He states that the purpose is 
to provide a continual monitoring of the amounts of wearmetal contaminants found in 
engine oil as an indicator of engine wear and possible engine failure [Ref. 2: p. 1]. The 
program technical manual defines the Joint Oil Analysis Program as: The Army, Navy, 
Air Force, Marine Corps coordinated effort to develop a standardized, mutually beneficial 
program to determine equipment wear condition through the use of oil analysis (Ref. 1: p. 
1-3]. 

The armed services were long under a preventive maintenance program which 
involved set maintenance at specific intervals. Although these programs were 
successful, they were costly both in terms of excessive time and costs of periodically 
replaced maintenance items. Within the world of lubricated engines, this meant an 
aircraft or vehicle might unnecessarily be taken out of service for maintenance. In 
response to this, on-condition or reliability centered maintenance was developed. At 
the heart of this maintenance idea is the performance of maintenance tasks as required, 
vice the arbitrarily set interval schedule. The end result 1s increased time between 
maintenance actions. Senholzi [Ref. 3: p. 18] addressed the benefits of increased 
intervals to include reduced maintenance costs, reduced equipment downtime, and 
reduced safety hazards that may occur between maintenance intervals. The nature of 
the oil analysis program ideally suits itself to this modern approach to maintenance. 

The oil analysis program uses used oil samples for testing and thereby is able to 


non-destructively monitor the apparent contaminant content in a lubricating system. 


This is done by measuring the wearmetal concentrations of certain key elements within 
an oil sample. The current values of this wearmetal, and the recent trend of those 
values are important in ascertaining the possible condition of a system. Through the 
use of the applicable technical manual, high levels and problem trends may be 
discovered. Determination of what action is to be taken and areas to check or 
continue to monitor are recommended by the technical manual. Unfortunately, the oil 
analvsis program is not a cure-all and will not discover every type of failure or 
guarantee to find a problem trend every time. 

The oil analysis program has vastly increased over the years as economics and 
costly, sensitive components have played larger roles. Currently, the primary oil 
analysis instrument is the A/E35U-3 (hereafter called the Dash-3). It is in use 
worldwide by the armed services and other groups, and plays a large role in the 
maintenance previously described. This is demonstrated by the 1.5 million samples 
analyzed per year by the 100 plus Air Force JOAP laboratories alone [Ref. 4]. Due to 
its weight, size, and supporting peripherals required, the Dash-3 is normally used in 
stationary laboratories (including aircraft carriers). 

In many areas, and for many engine types, the need for analysis is immediate in 
order to keep planes in the air and vehicles on the road. Restrictions on program 
usefulness are therefore evident if oil samples were required to be shipped to distant 
laboratories for analysis and have the corresponding delay in return of the subsequent 
results. The resulting requirements for a transportable oil analysis machine led to the 
Portable Wear Metal Analyzer (hereafter called the PWMA). It is a lightweight 
instrument that is easily transported utilizing two small cases. Its design makes it ideal 
for units that are deployed away from a nearby JOAP laboratory and which still need 
quick, reliable results. In the following sections, the principles of operation of both 
instruments will be discussed. 


B. DISCUSSION OF INSTRUMENTS USED IN ANALYSIS 

The measurements of two different types of spectrometers are compared in this 
analysis. The difference in their measurements is rooted in the different physical 
principals they employ, atomic emission and atomic absorption spectrometry. The 
following is provided as background for the physical concepts and uses of these two 


instruments. 


Pz 


1. A/E35U-3 

The A/E35U-3 (Dash-3) fluid analysis spectrometer is the standard instrument 
for the Joint Oil Analysis Program today. It is an atomic emission spectrometer and is 
manufactured by the Baird Corporation of Bedford, Massachusetts. It is a self- 
contained unit that weighs 800 pounds and has dimensions of 60 x 52.5 x 49.5 inches. 
Within its shipping container, it weighs 1400 pounds and has dimensions of 74 x 67 x 
47.73 inches [Ref. 1: p.4-2]. The instrument has the capability of measuring twenty 
elements within an oil sample. These twenty elements are listed in Table 1 . The 


measurement of all elements is accomplished simultaneously. 


TABLE 1 
PPE MEN S@fEAS@RED BY THE A/E35U-3 


Fe Iron 2 Sodium Ba Barium 

Ag Silver Ni Nickel Cd Cadmium 

AI Aluminum Pb Lead Mn Manganese 
Be Bervllium Si Silicon Mo Mel: denum 
Cr Chromium Sn Tin | V Vanadium 
Cu Copper Ti Titanium Zane 

Mg Magnesium B- Boron 


Atomic emission spectrometers are optical type instruments used to determine 
the concentration of wearmetals in lubricating fluid. The analysis is accomplished by 
subjecting (“burning”) the sample to a high voltage spark which energizes the atomic 
structure of the metallic elements, causing the emission of light. Within the Dash-3, 
this is done through the use of a rotating disk electrode spark source. The light is 
caused by the atomic structures of the elements being excited into higher energy states 
by the high temperature and then returning to their ground state through the emission 
of light. One or more characteristic wavelengths may be given off by each element. 
The emitted light is subsequently focused into the optical path of the spectrometer and 
separated according to wavelength, converted to electrical energy and then measured. 
The emitted light for any element is proportional to the concentration of wearmetal 
suspended in the lubricating fluid. This makes possible the estimates for the wearmetal 
concentration levels for each element that are normally measured in parts per million 
(ppm). [Ref. 1,5: pp. 2-1, 21] 

The final values obtained from a spectrometer are subject to many sources of 
error. Several of these sources apply to both the Dash-3 and the PWMA. Some of the 
reasons for error and factors that may influence the actual value obtained are given 
below [Ref. 1: pp. 2-3,2-4]: 


NS 


¢ Oil sample from wrong component 

® Contaminating substances in system 

e Oil is not homogeneous in content 

© «Type of instrument being used, emission or absorption spectrometer 

¢ Contaminated/poor quality calibration standards 

e Changes of rod or disc electrodes without restandardizing equipment 

e Metallic components in fluid additives within oil manufacture 

e §=Electrolytic corrosion due to seasons, storage 

e Dirt and sand contamination 

e Metal content of the fuel 

e Break-in period of new or recently overhauled components 

Another possible problem source is that of using synthetic ester oil to 
calibrate the instrument. Rhine [Ref. 6: p. 39] found this to directly effect the Dash-3 
by the enhancement of emission from metals in an ester oil matrix relative to metals in a 
hydrocarbon matrix. The end product of this 1s superficially high readings from the 
ester oil when calibration is done with a hydrocarbon oil. Similarly, problems are 
introduced by calibrating an instrument using standards prepared in a synthetic oil 
matrix. Specific to the Dash-3 design, errors may be introduced if there are 
inefficiencies in the particle transport capability of the rotating disk electrode, or in the 
vaporization and excitation of those transported particles at the arc/spark source 
(Ref. 6: p. 53}. A study for the Air Force by United Technologies Corporation placed 
some of the precision problems of the Dash-3 on the variable excitation characteristics 
of direct current (dc) arcs operated in air, inaccuracies due to variations in sample 
viscosity, and the inability to analyze particles larger than 10 micrometers in size 
[Reme7: p.1). 
2. PWMA 
The Portable Wear Metal Analyzer (PWMA) is planned as the standard 

instrument for deployments to remote areas where the Dash-3 is not located or 
suitable. It is a graphite furnace atomic absorption spectrometer manufactured by the 
Perkin-Elmer Corporation’s Applied Science Division of Pomona, California. It is 
composed of two self-contained units that double as shipping containers. The units are 
lightweight, weighing 60 pounds and 40 pounds respectively, and have identical 
dimensions of 11 x 18 x 18 inches [Ref. 8: p. 1-4]. The instrument has the capability of 


measuring nine elements within an oil sample, these elements are listed in Table 2 . 
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The measurements of all elements is accomplished simultaneously. The instrument has 


repeatability specifications of +/- lppm or +/- 2.5% of full scale, whichever is greater. 


TAGES? 
ELEMENTS MEASURED BY THE PWMA 


Fe Iron Cr Chromium Ni Nickel 
Ag Silver. Cul Copper Si Silicon 
AT Aluminum Mg Magnesium Ti Titanium 


Atomic absorption spectrometers are similar to their emission counterparts in 
being optical type instruments used to determine the concentration of wearmetals in a 
lubricating fluid. However, the physical concepts of the machines are very different. 
With the atomic absorption spectrometer, the atomic structures of the elements present 
become energized enough by high temperatures to absorb light (rather than emit light 
as in the Dash-3). Within the PWMA, the heating is done in a graphite furnace tube 
whereas other atomic absorption spectrometers usually use a flame. The graphite 
furnace atomic absorption spectrometer uses the graphite tube as a resistive heating 
element to atomize the wearmetals. Light energy having the same (“characteristic”) 
wavelength of the element being analyzed is then radiated through the graphite tube. 
Two multi-element hollow-cathode lamps are used for this in the PWMA. The 
resultant light is converted to electrical energy and measured electronically by using 
photomultiplier detectors. The amount of light energy absorbed by the elements in the 
optical path of the graphite tube is proportional to the concentration of wearmetal 
suspended in the lubricating fluid. This makes possible the estimates for the wearmetal 
concentration levels for each element that are normally measured in parts per mullion. 
The after to before ratio of light intensities that is used makes the signal less 
susceptible to intensity variations of the spectra caused by the environment. 
[Retee4: pp 2-2,6] 

Six prototype PWMA instruments were acquired by the Air Force for testing 
in an operational environment with military personnel as operators [Ref. 9: p. 2]. Field 
test plans were developed by the Air Force Wright Aeronautical Laboratories 
(AFWAL). As part of this plan, instruments were provided to four locations to be 
used alongside the Dash-3. Locations in receipt of these instruments were: Naval Air 
Rework Facility (NARF), Pensacola; Elmendorf AFB; Langley AFB; and Myrtle Beach 
AFB. Preliminary testing on a prototype PWMA showed it to be successful in 


ls 


measuring particles up to 20 micrometers in size [Ref. 4: p. 60]. The prototype 
PWMA’s currently use preprogrammed algorithms in calibrating themselves. This is 
done using three successive calibration runs. The runs specify three points for each 
element’s calibration curve algorithm and may result in small inaccuracies. Further 
areas of possible error are described in Reference 4. 

Initial testing of the PWMA has shown it to yield values that are 
approximately one-half of those obtained with the Dash-3. For example, a sample that 
measures 11 ppm on the Dash-3 may only measure 6 ppm on the PWMA. The 
wearmetal concentration level (in ppm) is readily seen from this to depend on which 
instrument is used. Investigation of this point shows the controlling factor to be in 
whether the measurement is based on atomic emission (as with the Dash-3), or based 
on atomic absorption (as with the PWMA). This difference in using light emission or 
light absorption for estimating an element concentration level leads to the separate 


scales obtained. 
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Il. DATA 


A. FIELD USE OF JOAP DATA 
Once values have been obtained from an oil analysis spectrometer, an analyst 

must be able to do more than add them to a database. His primary goal is to 
determine if there is any abnormal condition present or developing in the parent 
component of the sample. This goal may be achieved by following the wearmetal 
methodology presented by the JOAP Laboratory Manual. The major parts of the 
method are {Ref. 1: p.6-1]: 

¢ awearmetal evaluation and trend table 

e atable providing decision making guidance 


¢ supplemental information providing additional wearmetal diagnostic guidance for 
each type of equipment evaluated 


When consumption and replenishment of fluid is continuously balanced in a 
component, the wWearmetal concentration within the lubricant tends to increase after a 
fluid change until a steady state level is achieved. This effect is shown in Figure 3.1 . 
The wearmetal concentration steady state point in this model is thus a function of two 
variables: rate of fluid consumption and replenishment; and the rate of wearmetal 
production by internal friction within the equipment. When the fluid is added or 
changed periodically, as is done with many items of equipment, a sawtooth pattern 
emerges. Such a pattern is shown in Figure 3.2. {Ref. 1: pp. 2-5,6] 

Because of a tendency toward a steady state level, it is not only important to 
know the actual wearmetal concentration, but also to know the rate of increase (trend) 
of that concentration. Increased sampling is always recommended as levels increase in 
order to decrease the probability of a wearmetal concentration going from normal to 
abnormal without an intermediate sample. The exact procedure the oil analysts are 
trained to follow in evaluating a sample is a logical sequence of actions. The JOAP 
manual mandates the following procedure [Ref. 1: p. 6-2]: 

e §=©Analyze the sample and obtain the wearmetal concentration results. 


e Determine the range where each critical wearmetal falls from the appropriate 
wearmetal evaluation criteria and trend table. 


¢ Compare the wearmetal concentration ranges of the current sample with the 
ranges of the last previous sample. 


e Determine the wearmetal trend (trend value for a ten hour period). 
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e lf a maintenance action ts indicated, consult the supplemental information 
provided. 





Figure 3.1 Theoretical Concentration with Constant fluid Replenishment. 





HOURS 
Figure 3.2 Concentration with Periodic Fluid Addition or Change’. 


Table 3 [Ref. 1: pp. 98,99] is provided to gatn a fuller understanding of the oil 
analysts’ procedure. ft provides the wearmetal evaluation criteria and trend table for 
the 164 series engine commonly used in the CEE-53A,D helicopter. A sample scenario 
would show the analyst checking for a normal trend by comparing the current values 


with those of the previous sample, and then checking for to sce if the current values are 
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in a normal range. With that information he would-enter a decision making guidance 
table and obtain from it the advice codes it lists. This decision aid is shown in Table 4. 
When the code indicates anvthing abnormal, some type of action is recommended. 
Normally the initial action is to request a ‘red-tagged’ sample and to cease operation of 
the equipment. The ‘red-tagged’ sample is to gain confirmation of the abnormal results 
before any tvpe of maintenance action will be recommended. If abnormal results are 
confirmed, the supplemental information provided in the technical manual is consulted 
and a maintenance action is recommended. The supplemental information may be very 
useful in pin-pointing possible problem areas and identifying a specific recommendation 
to the sample originator [Ref. 1: p. 6-1]. 

An example of the above procedure using the tables shown is given here. An oil 
analyst receives a sample from an engine that was previously sampled twenty engine 
hours earlier. His logbook shows that its last sample was normal and had values as 
follows (in ppm): Fe-10, Ag-0, Al-2, Cr-0, Cu-1l, and Mg-l. The present sample is 
tested with the following results: Fe-19, Ag-4, Al-6, Cr-0, Cu-1, and Mg-l. The analyst 
determines the ten-hour trend values to be (in ppm/hour): Fe-4.5, Ag-2, Al-2, and the 
remainder zero. Using Table 3, he would discover Fe and Ag to have abnormal trends 
and high ranges, with all other elements normal. Entering Table 4 in the high range in 
column one and reading across through the normal range in column two for the 
previous sample, and then through the abnormal trend in column three, one obtains 
the advice code P from column four. The analyst translates this code as: do not 
fly/operate; do not change oil, submit sample ASAP and relays this to the sample 
Originator. The ‘red-tagged’ sample is subsequently analyzed for verification of the last 
results. If results are identical to the previous ones, the appropriate advice code from 
column five is utilized. This might include using the lower part of Table 3 to determine 


possible maintenance problems such as the power turbine shafts. 


B. DISCUSSION OF DATA USED 

The data in this analysis was obtained from a field test for the evaluation of the 
Portable Wear Metal Analyzer that was conducted in 1985. The data was collected 
over a two month period from four locations: NARF, Pensacola, Elmendorf AFB, 
Langley AFB, and Myrtle Beach AFB. Both Dash-3’s and prototype PWMA’s were in 
the laboratories at each site in order to provide data for their comparisons. Only the 
nine elements that the PWMA is able to measure were recorded for the Dash-3. 
Though the primary objective of the field test was to determine the functionality of the 
PWMaA [Ref. 10: p.1], the data obtained easily lent itself to this study. 
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Ten ts3 
SAMPLE WEARMETAL EVALUATION CRITERIA AND TREND TABLE 


ENGINE: T64-6/-7/-7A/-10/-12/413/415 
AIRCRAFT: (CH-53A/D) (RH-53A/D) (VH-53D) (UC-8A) 


A/F.35U-3 
Abnourmai Trend 


(PPM Increase 
in 10 hes) 


Marginal Range 15-17 


Abnormal 





Average Concentration Other Elements: 
Niel Po=2 Si-3 Sn°7 Tiel Mo-1 


ENGINE: T64-6/-7/-TA/-12/-4 13/415 (Cont) 
AIRCRAFT: (CH53 A/D) (RH-53A/D) (VH-53D) 


Main bearing bails, roliers and races, accessory gearbox gears and shafts, PTO gears 
Ag Accessory gearbox bearings 
Power turbine sha(ts 
PTO beanngs 
A Cu Main bearing caves 
Lube and scavenge oil pumps and thermal gradient housing 
NOTE 


Thermal gradient housing could be pnme source of [Al] especially after heavy use of 
engine anti-icing system. 


TABLE 4 
DEGISTON MUAKING GUIDANCE TABL] 


~ 
4 
A 
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As designed, the test called for each site to use six graphite tubes for 160 sample 
‘burns’ each. Many samples were to be tested five consecutive times when analyzed. 
The test plan was designed for repetitive four day cycles. The sequence of events it 
followed for the PWMA is shown in Table 5 [Ref. 10: p. Il]. The sequence of events 
for the Dash-3 was nearly identical to that of the PWMA. Consistency in following 
this test plan varied among the different locations. This directly affected the number of 
sample results provided by each location. A compilation of the number of samples 
used in this analysis from each location is given in Table 6. 


TABLE 5 
PWMA FIELD TBs sBURN SEQUENCE 


Day | Calibration Runs 
Verification Samples 
orrelation Samples |, 2, 3 
Random Samples 
Verification Samples 


Day 2 Reslope | 
Verification Samples 
Correlation Samples 4,5 
Random ae 
Calibration Check Samples 
Random, Samples 
Verification Samples 


Day 3 Calibration Samples 
Verification Samples 
Correlation Samples 6,7 
Random, Samples 
Verification Samples 

Day 4 Reslope . 
Verification Samples 
Correlation Samples 8,9 


Random Sa 
Calibration Check Samples 


As noted in Table 5, there were four main classes of samples tested with the 
instruments: calibration, verification, correlation, and random samples. Calibration 
samples of MIL-L-7808 oil were used for the internal calibration of the PWMA. Three 
calibration samples containing 20%, 50%, and 100% respectively of full scale values 
for all nine metals were used for the calibration runs and were also used for quality 
assurance in the calibration check runs. Verification samples were made in MIL- 
L-7808 oil and contain concentration levels of 10%, 40%, and 70% of the top of the 
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TABLE 6 
NUMBER OF SAMPLES USED IN ANALYSIS 


Elm Lan Myr Pen 


Correlation Samples 3 26 ee ae 
| Verification Samples 65 Li, 39 


dvnamic range. Correlation samples were designed to reflect the pseudo-realism of 
having various concentration levels for each element. Nine different correlation 
samples were used in the test. Because of the inavailability of MIL-L-7808 oil with 
significant wearmetal levels for the correlation samples, used MIL-L-23699 turbine 
engine oil was utilized after blending with organo-metallic concentrates as necessary by 
JOAP-TSC. These samples were ensured to contain a minimum of seven wearmetals. 
The random samples were routine daily samples obtained at each site and were 
typically MIL-L-7808 lubricants from the Air Force, and MIL-L-23699 lubricants from 
the Navy. Reslopes were calibration checks using a 50% of full scale calibration 
sample. [Ref. 10: p. 6] 

Due to the nature of the calibration runs and the low values in the random 
samples, only verification and correlation samples were used in this analysis. This 
made possible a total of twelve samples, consisting of three different verification and 
nine different correlation samples. Sample numbers (or levels as appropriate) were 
kept in the database as well as the graphite tube number in which a sample was used 
for later use in determining their means. Each sample ‘burn’ with the PWMA 
correlation samples was repeated five consecutive times. In order to provide a one to 
one correspondence of data between the instruments, the average of the consecutive 
‘burns’ was used in the analysis. The PWMA value (average value for correlation 
samples or actual measurement for verification samples as appropriate) was coupled 
with its counterpart on the Dash-3 that was closest in time. In normal cases this led to 
results such as pairing the first verification sample ‘level 10° measured during the day 
from each instrument. Differences in sample sizes for the various element-location 
combinations was a result of missing data points and the removal of extreme outliers 


from the data. 
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IV. ANALYSIS 


A. OBJECTIVES 

The objective of this analysis is to find a functional relationship between Dash-3 
and PWMA measurements for the nine elements that the PWMA measures. Analysis 
of variance techniques will be used in order to determine what effects cause significant 
differences in the measurements obtained. Regression analysis will be used to yield the 
parameters of an equation to relate the instrument measurements, and to obtain the 
standard deviations of those parameters. The coefficient of determination or r? statistic 
will be used as a measure of how well the model fits the data and to compare the 
different models. It is also Known as the square of the coefficient of determination and 


represents the percentage of total variation explained by the model. 


B. SCATTERPLOTS 

Scatterplots of the data are useful in determining areas for further investigation. 
They were used in this analysis to provide insight into the rejection of the hypothesis 
that there is no difference 1n measurements due to location. The rejections are 
described in the analysis of variance section. The scatterplots in Figures 4.1 - 4.9 are 
location coded and include all verification and correlation sample data points. As 
reflected in Table 7, examination of these plots show Pensacola with higher values and 
Elmendorf with lower values for most elements (seven of nine in each case). The 
consistent extremes shown by these two locations are a primary reason for differences 
found between locations. The number of samples submitted does not explain the 
extremes as both had larger sample sizes. Coded scatterplots of the graphite tube 
furnaces for each location/element combination are shown in Appendix A. The tubes 
were reviewed in particular because they were prototypes and not production models. 
The quality controls that will be used in manufacturing the production models was not 
used for the prototypes, consequently, it was suspected that they could be a large 
source of variability in the measurements. The graphite furnace tubes were provided to 
each site with the tubes labeled sequentially within sites. Some tubes malfunctioned 
and could not be used due to calibration failures, causing skips in the tube numbers 
reported. Most tubes did not last as long as expected originally so all tubes did not 


have an equal number of burns. With minor exceptions, no unusual characteristics 
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Location Coded Scatterplot of Fe Using All Points. 
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Figure 4.2 
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Location Coded Scatterplot of Ag Using All Points. 
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Figure 4.3 Location Coded Scatterplot of Al Using All Points. 
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Figure 4.4 Location Coded Seatterplot of Cr Using All Points. 
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Figure 4.5 Location Coded Scatterplot of Cu Using All Points. 
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[igure 4.6 Location Coded Scatterplot of Mg Using All Points. 
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Figure 4.7 Location Coded Scatterplot of Ni Using All Points. 
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Figure 4.8 Location Coded Scatterplot of Si Using All Points. 
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Figure 4.9 Location Coded Scatterplot of Ti Using All Points. 


TABLE 7 
SUN ai OF EOCATION CODED SCATTERPLOTS 


Element Location Coded Results 


[e low 
Ag | P high 
Cr 
\tg P high 
Pe Mei ae 


were found. A summary of the graphite furnace tubes in Table 8 shows Elmendorf 
used six different graphite tubes in its experiments with tube I measuring higher than 
normal in seven elements, tube 7 high in five elements. and three points of tube 5 high 
in seven elements (indicated by 5,). Though the three points from tube 5 appeared to 


be an aberration, they were included in the analvsis. Higher,lower than normal ts not 
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TABLE 8 
SUMMARY OF GRAPHITE TUBE CODED SCATTERPLOTS 


a 

Te [asym [mine | 

a 

ee 

a 

i 
[tie 
















used in a statistical sense but as a subjective term. A symbol that characteristically 
plotted as a maximum/minimum value across the range of x values was given the 
appropriate subjective term. The five tubes used by Langley showed tube 7 measuring 
low in three elements, tube 2 low in one, and tubes |, 3, and 5 high in one element 
apiece. Myrtle Beach used six tubes and showed nothing unusual except for tube 3 
measuring low in one element. Pensacola provided the most data and used the most 
tubes, 11. Its scatterplots are not unusual except for tube 1 being high for two 
elements and low for one,,and tube 16 being high for one element. In summary, the 
scatterplots indicate that the locations provide different measurements, with Pensacola 
normally high, and Elmendorf normally low. The graphite furnace tubes do appear to 
have a Statistically significant difference within each instrument, but not a practical 


difference. 
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C. ANALYSIS OF VARIANCE ; 

As stated, a superior tool to scatterplots in determining if various effects are 
statistically significant is analysis of variance. It is used here to determine if the 
Various instruments (the four prototype PWMA’s and the four Dash-3’s) are providing 
measurements close enough to be considered not statistically different. Examples of 
the various ANOVA models applicable can be found in References 11 and 12. Using 
the standard assumptions of normality, an alpha level of significance of 0.05 is used 
throughout this paper. Statistical significance obtained may be more restrictive in its 
determinations than what would be considered a practical significance by the oil 
analysts. In other words, the numbers may be close enough to be considered the same 
by the analysts. No assertions as to what this practical significance might be will be 
addressed. Analyses of variance were run on the correlation samples to determine if 
location, sample, or their interaction is significant. The resulting tables and associated 
p-values are given in Appendix B. As summarized in Table 9, location differences 
were found with all elements except Ni on the Dash-3. Location-sample interactions 
were significant with four elements for the Dash-3, and in all but two elements for the 


PWMA. These results imply greater variability in the prototype PWMA instruments. 


TABLE 9 
SUMMARY OF CORRELATION DATA ANOVA 





i i X x x X 
ee <P Eee Sea Xx 


X denotes statistically significant at a= .05 


The analyses of variance on the verification sample data for both instruments 1s 
provided with their corresponding p-values in Appendix C, and a summary given in 
Table 10. Table 10 shows that when interactions were considered, location was not 


oil 


significant for three PWMA elements (Cu, Si, and Ti) and one Dash-3 element (Ni). 
Analysis of interactions showed all PWMA elements except Mg did not have 
significant interactions, while all Dash-3 elements but Ni did have a significant 


interaction. In contrast to the correlation data, this shows a strong variability in the 


TABLE 10 
SUMMARY OF VERIFICATION DATA ANOVA 


PWMA Dash-3 





X denotes statistically significant at a= .05 


Dash-3 compared to the prototype PWMA. As expected and designed, the different 
verification levels and different correlation samples were significantly different within 
their respective groups. Further analysis on the verification data is summarized in 
Table 11 , which takes a closer view of location effects at each of the three levels. 
Recall that these levels represent the samples which contain 10, 40, and 70 percent of 
the full scale dynamic range for each element. The table gives each instrument's 
corresponding element-level matrix of significant location effects. The matrix shows 
the prototype PWMA to be very consistent in its measurements across locations as 
Opposed to the Dash-3, which lacks consistency for most elements and levels. The 
basis for this statement lies in the few holes shown in the Dash-3 matrix and the few 
marks in the PWMA matrix. Examination of location effects on the verification data 
was also done using means. The highest and lowest location mean was recorded for 
each element and is shown in Table 12. From this it can be seen that a location that 
is high(low) for an element is usually high(low) for all three levels for that element. Of 
interest is that the location that measures high(low) for an element on one instrument 


will not normally be high(low) on the other instrument as well. This shows an 
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TABLE lL 
SUMMARY OF VERIFICATION LEVEL ANOVAS 





X denotes statistically significant at a= .05 


apparent low correlation between the two instruments regarding the same location 


being high(low) for a particular element. 


PARE a 
VERIFICATION LEVEL HIGHEST AND LOWEST MEANS 


PWMA Dash-3 


S1 i P P 
E ACEI 
P 
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D. REGRESSION 

Regression is used as the principal tool to determine functional relationships 
(transformation equations) between the measurements of the two different instruments. 
Tables 13 - 21 show for each element the results of several regressions run on the data. 
For each regression, the following are given: r’, the intercept and its standard 
deviation, the slope and its standard deviation, and the number of data points used in 
the regression. Two principal models were analyzed, the first model used all data 
points from all locations regardless of the location sample sizes (model 1); and the 
second model used the means of the samples and levels at each location, thereby 
allowing equal sample sizes (model 2). The various sample sizes were presented in 
Table 6. Appendix D gives a basic regression model. The majority of regressions 


done used the simple linear model of: 
PWMA = Bo + B,(DASH-3) + e. (eqn 4.1) 


Other regressions used a quadratic fit or used simple linear fits on transformations of 
the data. No attempts were made to force the resulting equations through the origin. 
An ideal regression using the simple linear model with least squares estimators 
has a high r* . It will also meet the assumption of constant variance across the range of 
interest (homoscedascity), and that of normalized residuals. The usual least squares 
regression estimators have optimal properties if one assumes that the variance of the y 


values remains constant as X increases, together with 
Ely. a> 2 sxe. (eqn 4.2) 


Many measuring instruments have variances which increase with the magnitude of the 
quantity measured. Letting 6? represent the variance of y,, one could assume a simple 


linear relation 
Co KX, , (eqn 4.3) 


in which case Equation 4.2 can be transformed by dividing by JX, . This will give a 
constant variance model and help to normalize the residuals. However, this procedure 
will not work if the data displays a variance that first increases and then decreases as 


the true content increases. 
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TABLE 13 
REGRESSION RESULTS FOR FE 


Data 


oe 
Pesci | ruin | os [as [as [oor | ow | ni 
Pattee rin at fae [ae Po [oa 
Tonite [ite [oe [ae [ow [or fas [ar 

ae = | 


“Geves [aoe oe pie Par fom | 
Ce EC 
Seen — [oe oe [ae Pin [aD De 
pia [ait a | = [= >? | - | 
Crna? [aes |= pp | 
puiatoans [att 7 |= [=f | |= 


*.denotes not included due to little or no r* improvement. 
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TABLE 14 __ 
REGRESSION RESULTS FOR AG 


Elmendorf | l- Lin | 


Coats [rs af Lae of 
“ineae [tum | ae | a | [ss [om | st 
Tarte [tin [| ae [as [si oe [ae 
Cisne [tin ar [ae as [ons | 
raniee Peed ae [Ds ae 


even (ia a fo [oe foe 
Teatvieae [atin [ os [as | ae | oo fom [a 
Catterwen [oun | sane [so [os [a 
Prwwiatoee [aim[ f= p= >= [+ [a 
Crwwapans® [aim| asf p= p> | | @_ 

pwaatbans [aul |= [= ]- |= [a 


*.denotes not included due to little or no r* improvement. 
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TABLE 15 
REGRESSION RESULTS FOR AL 





Elmendorf l- } i-Lin 


Se ee 
sont [ie [ ow fae |e [aw [os 
Crenacse pun | 9 fs | vo fom [on [us 
Pattee [run] or [as Dw [ae [oor | 
tannin [tin | ai [7a] ss | as | mas | a 
a 


ceria pee) @ a pw pes pes pa 
veatviene [2m | ow [100 [so | mp om [ae 
puNabane [atm | ~ |= | | | «_ 
Trwwadans® [aim[ of)? [- | «_ 
Prwwatowis [aim] [> [-]* [7 1«_ 


*.denotes not included due to little or no r- improvement. 
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TABLE 16 
REGRESSION RESULTS FOR CR 


Cissy fin | or [we | a | es [os os 
“sneer [in| 2 | | as [oe [oo lo 
essen [in| 96 | | ae | woe | 8 


rN 
ratios [reel se [TT 


ewes [simp oa | [aw fom |e 
[seater | aun] os [a | os [am | a | 2 
Taree wens [>in] 90 [06 | ae | me [ae | 
ce 
Cewwatoens [aunt of > [-]* [-]*_ 


*.denotes not included due to little or no r’ improvement. 
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TABLE 17) 
REGRESSION RESULTS FOR CU 





ease re 
anes in| eam [oe [oss | on | 7 
reco | Min] 90 | 266 | | 7 [om |e 
ratte | vin | [ais [a6 [or [ow [a0 
iawyres | in| [an | @ [7a [on | a0 
partes [eed ae [P= Te 


ceeves [am | a ae je [a pap 3 
[veafviens [atin] [an | «| oe [ow | a 
Taitoe venetian [ oe [1s [vo [ose [om | as 
Tinie Mean [atin | oa [aw [re | [om | a6 
rwviatpaaa? faumf af 1 -[- | |*_ 
Trwwabans [aiml f= ]-]- [+s 
Prwwatoana [auml f= [-{- |-1s_ 


*.denotes not included due to little or no r* improvement. 
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TABLE 18 
REGRESSION RESULTS FOR MG 


| Elmendorf 1- | l-Lin ae 
a 
Montedack [utin| #2 [207 [us | se | oo | ot 
Tatice | vtin[ te [aw | © | se [oe | 29 
Ttanvirren | ein | as [as | o> | our | om | 
Pantee [eed = Te To 


reeves ppt) we [aa pi aw pow | 
(vest vieane [aim | so [aes [aie [os [on | a 
Tairtoe wes [aim |e [206i [ser oo | os 
Cisir wean [aun | a [usr [iar [on | ows | a4 
Prwwiatoeme [pim| |= p=] * [+] «_ 
Crwwabans® [aim |] > p= [> | | 
Prwwatpeas [amt p= f= f= 1-1«_ 


*.denotes not included due to little or no r improvement. 
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TABLE 19 
REGRESSION RESULTS FOR NI 


| Elmendorf - SI-Lin | 80 | 


“ts Tia a ae Dao D9 
“Menieoach [tun | os [aes [st [a | me | 0 
Pec pase oe ae [a a Do 
Nite | in| [| a0 [| 0 | oe 
conten [vain [os ae [st [om 
rane feed 


—— Means pis at Te fe fom Ta 
VerifMeans | 2tin | 98 | 3.68 | 70 | 793 | 032 | 12 
Pattee [ria] ow |» Lo [oe Po 
LEMP Means | 2tin| 97 | 37 | 42 | 638 | ois | 36 
Pwoa*Das-# | 2tin{ 39 | + | + | os | o* | 
| PwMA-Dest-s? | 2tin| a2 j * | * | * | + | 
jPwoa®Dashs | atin] 8s | + | + | * | * | 


*_denotes not included due to little or no r* improvement. 


1. Regression Using All Data Points Individually (Model 1) 
Seven different regressions were run where all data points from all locations 


were used. The first four of these were run by using a simple linear model on each 
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TABLE 20. 
REGRESSION RESULTS FOR SI 


[Biaeoe c= abet 

Caneny [um | ae [ss Pie Lae [oo 
“Mone bac [tin | ac | om [ ss [a7 fo | so 
Crenaw [rtm] a | 7 | @ | oo | me [us 
Tarte [vim] [ar | [as [on | a 
Tianna | rtim| oo [ae | ae [ane [ow | 
rae [ve Pe | 


Coven pate] @ [a [= [oe pe pa 
eaters [oun | o« [un [ia [0 | oo | a 
Tatton [ain as [ot | ef [on Pa 
ive [aim [a7 [0 | [ae om | a6 
Crwiatomn [pam | art = p>) > [a 
Trwwaban® [aim] |? [= [> [= | @” 
Trwwatduna [aim] @ | [>> | 1 


*.denotes not included due to little or no r’ improvement. 

















location and the data then used to gain a parametric understanding of the scatterplots 
presented earlier. Non-parametric results in Tables 22 - 24 were derived from these 
regressions. Table 22 ranks the slopes for each element by location and shows 
Elmendorf to have a smaller slope with six of the nine elements (explaining the earlier 
scatterplot results). It also shows Pensacola to have a higher slope than the other 
locations in five of the elements. The rankings of the standard deviations of the slopes 
are given in Table 23. It was found to be smallest most often with Pensacola, with 
Elmendorf exhibiting the second smallest standard deviation. Table 24 ranks r? and 


2 


shows Elmendorf to have a lower r“ in seven of the nine elements. 
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TABLE 21 
REGRESSION RESULTS FOR TI 


| Elmendorf | A-Lin | Elmendorf | A-Lin Lin a 

Langley “ap ae ae 
Myrtle Beach sttin| 92 | 11s | 4a | 648 | 027 | 50 | 
Press [eae DD or [a | 
Altec | ttn | si | sar | 25 | 602 | ore | 343 
Cease | rin) ao) aoe fos | a 
pas [re = 


Corr Means ee 
ee | Oe 
& All Loc Means | 2Lin | 91 | -.42 4 | 607 | 028 | 48 | 
Sister [ae [oe [oD fo or [ 
| PWMA?-Dash-3? | 2-Lin Catlett] 8 
Crennosest [ae f= P= T= Ta 
Trwwatdens [atin] [> |]. 11s 


*.denotes not included due to little or no r? improvement. 


The remaining three regressions were used to determine the best way to fit the 
data from the combination of locations. A combination using all locations except 
Elmendorf was investigated because of Elmendorf’s significantly smaller slope. For 
most elements, (all except Al and Mg), increased when Elmendorf was not included. 
No more than a .02 improvement was gained by using a quadratic regression on the 
data. Because of this small gain, it was decided not to pursue this regression model. 
Other regression variations were tried in the sample and level means model described 


next. 
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TABLE 22. . | 
RANK OF SLOPES FOR INDIVIDUAL LOCATION REGRESSIONS 





Rankings highest(1) to lowest(4) 


TABLE 23 
RANK OF SLOPE STANDARD DEVIATIONS FOR LOCATIONS 





Rankings smallest(1) to largest(4) 


TABLE 24 
RANK OF R* FOR LOCATION REGRESSIONS 


Tal alae. aps] a [a 
Pele Te fe 
habbit. 
reli bpp 
reef te pe} pe pepe [sf 


Rankings largest(1) to smallest(4) 
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2. Regression Using Sample and Level Means (Model 2) 

Due to the disparity in the number of samples provided by each location, 
(from a low of 51 to a high of 119), a method was needed to eliminate the skewing of 
results toward the locations with more data points. By taking the means for each of 
the three verification levels and of the nine correlation samples, all four locations 
would be equally weighted. This approach allowed a total of twelve means (9 
correlation and 3 verification) to be provided from each location. Several regressions 
Were conducted using this approach. Linear regressions were conducted on the 
verification level means alone, the correlation sample means, on the combined means, 
and on all means except those from Elmendorf. Linear regressions on transformations 
of the combined location model were also investigated. From Tables 13 - 21 , it can be 
seen that the slopes for verification level regressions were often high in comparison to 
the other regressions. The low number of data points (three clusters of four points 
each) probably contributed to this. The correlation sample means regression was non- 
descript and quite similar to the regression utilizing model | for many of the elements. 
The model 2 regressions using the mean values from all locations and using mean 
values from all locations but Elmendorf are also very simular to their counterpart model 
1 regressions that utilize all points. The similarity can be seen in the summary of 
equations table presented in Table 25. The equations shown have the highest r? for 
each respective model. The supporting data either consisted of data from all locations 
or all locations except Elmendorf (indicated by LMP). The r values using the means 
(model 2) is greater than or equal to those obtained using all points (model 1) in every 
instance. Improvements to the r values were investigated through three linear models 


which used transformations of the data, these were: 


PWMA? = By + B,Dash-3? + e, (eqn 4.4) 
PWMA = Bo + B,Dash-3* + e, (eqn 4.5) 
PWMA? = Bo + B,Dash-3 + e. (eqn 4.6) 
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Throughout all of the regressions, a recurring problem of violating basic model 
assumptions was observed. This is not to imply the equations are not unbiased 
estimators, but only that they may not be the best estimators. 

The common data transformation technique described at the beginning of the 
section was attempted but it had little success in normalizng residuals or gaining 
constant variance. Four elements (Al, Cr, Ni, and Ti) did attain a constant variance, 
and five elements (Fe, Al, Cr, Mg, and Ti) achieved normalized residuals. Only the 
elements Al and Cr met both assumptions of constant variance and normalized 
residuals. With only two elements meeting both assumptions and four others meeting 
only one after using the transformation derived from Equation 4.3, it was decided not 
to pursue the transformation farther. This decision was supported by taking into 
account the small improvement gained in meeting the assumptions at the expense of 


the complexities that would be introduced into the equations for the oil analyst. 


Wea 2) 
REGRESSION EQUATIONS FOR MODEL | AND MODEL 2 


oni Tes [Toes Tae 
Te pena am | Ewe) a [| Pa-2004 6D 

se [Po lon 090 roar arr 308 ow 
Car [p= sisrae> [als |[ po oreramad | LMP] 
Cee [r= aera | owe 9 || pon aeremp | ume] ot 
rer | pearsenap | owr[ a | po-aas.o9p | MP] oe 
Re ee 
CS [r-osrasi [ows [f= ave | uw] ow 
rs [pee aeasep | ow a [[ pos sorase | uw] 7 
ci [re ame [owe [fo osrs | uve] 95 
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V. SUMMARY AND CONCLUSIONS 


The Joint Oil Analysis Program has reached an important milestone with the 
addition of the Portable Wear Metal Analyzer. Strong relationships were shown to 
exist bv this analysis between the measurements taken by the prototype PWMA and 
the A/E35U-3. The analysis was based on data from a 1985 field evaluation of 
prototvpe PWMAs. Although the data was not obtained in a consistent fashion by the 
laboratories, it was useful in determining accuracy within and between the laboratories. 
A subsequent analysis similar to this one is merited when production PWMAs and 
supporting parts are in use. The equations below describe the best estimate of the 
relationship between the Dash-3 and prototype PWMA readings. The equations could 
be used in making conversions from the JOAP Laboratory Manual. Any values that 


result in a negative translation should be regarded as zero. 


TABLE 26 
RECOMMENDED EQUATIONS FOR A/E35U-3 TO PWMA CONVERSIONS 










PWMA = -2.09 + .617(Dash-3) 


ca [ria = 3+ a 
PWMA = - .29 + .822(Dash-3) 
ce R= a mo 


PWMA =_ .37 + .638(Dash-3) 
PWMA = -.80 + .356(Dash-3) 


PWIA = 2103 2 6la(Dash-3) 








The equations presented are unbiased estimators when a _ simple linear 
relationship is assumed (although they do contain violations of homoscedascity and 
normalized residuals). The ability of the equations to describe the data varied, but 


47 


seven of the nine elements had r* values above .80. (Al at .79 and Si at .77 were the 
exceptions). These values were obtained using prototype instruments that Were using 
prototype graphite furnace tubes and were therefore considered quite good. One of the 
four locations consistently had a slope less than the others, causing it to be dropped 
from the final analysis for seven of the elements. This improved the r* values in each 
of those elements. 

Items that may have influenced the data used included the non-production 
graphite furnace tubes, the particle size capabilities of each instrument, the nature of 
the correlation samples, and the several factors mentioned in Chapter II. The furnace 
tubes were found to have statistically significant differences at each of the locations, 
but they did not appear to have a practical difference for the oil analyst. Further 
analysis using production model tubes will serve to clarify this area. In regards to 
particle size, the Dash-3 can only measure particles less than ten micrometers in size, 
while the PWMA can measure particles up to twenty micrometers. It can be deduced 
from this that the PWMA/Dash-3 functional relationship is directly dependent on 
particle size. Therefore a sample with all particles less than ten micrometers in size 
should show a very strong relationship between instruments. However, it is possible 
for a sample to contain a majority of its particles in the ten to twenty micrometer 
range (or higher) and thereby yield readings not consistent with the equations derived. 
The correlation samples were manufactured using a blend of oils and enriched with 
powdered organo-metallic concentrates to obtain measureable wearmetal levels. The 
room for inconsistencies due to the type oil used for instrument calibration with the 
type in the correlation sample, and due to particle sizes in the samples leaves this as an 


area for further scrutiny. 
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APPENDIX A 
GRAPHITE TUBE CODED SCATTERPLOTS 


OF VERIFICATION AND CORRELATION SAMPLES 
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Figure A.l Graphite Tube Scatterplot-Location:Elmendorf-Element:Fe. 
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Figure A.2) Graphite Tube Scatterplot-Location:E]mendorf-Element:Ag. 
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figure A.3> Graphite Tube Scatterplot-Location:EImendorf-Element:Al. 
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Figure A.4 Graphite Tube Scatterplot-Location:Elmendorf-Element:Cr. 


CU IN PPM 
oS 
cO 
5 
b 1 *TUBE 1 
2 -TUBE 2 
Oo 
oO 3 TUBE 3 
' 4 -TUBE 4 
5 5 -UBE 5 
{ ' 4 
S Se 2 1 4 a ? =TUBE 7 
Qa. 





DASH- 3 


figure A.S Graphite Tube Scatterplot-Location:Elmendorf-Element:Cu. 
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Figure A.6 Graphite Tube Scatterplot-Location:Elmendorf-Element: Mg. 
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Figure A.7) Graphite Tube Scatterplot-Location:Elmendorf-Element:Ni. 
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Figure A.8 Graphite Tube Scatterplot-Location:Elmendorf-Element:S1. 
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figure A.9 Graphite Tube Scatterplot-Location:Elmendorf-Element:Ti. 
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Figure A.10 Graphite Tube Scatterplot-Location:Langley-Element:Fe. 
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Figure A.ll) Graphite Tube Scatterplot-Location: Langley-Element:Ag. 
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Figure A.!2.) Graphite Tube Scatterplot-Location:Langley-Element:Al. 
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Figure A.13) Graphite Tube Scatterplot-Location:_Langley-Element:Cr. 
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Figure A.14 Graphite Tube Scatterplot-Location:Langlev-Element:Cu. 
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Figure A.IS Graphite Tube Scatterplot-Location: Langley-Element: Mg. 


56 


NI IN PPM 


> elUAF 2? 


50 


3 °fUDE 3 
5 +1UBE 5 
1 


y -fUGE 7 


FWMA 





0 20 40 60 
bASH—3 


Figure A.16 Graphite Tube Scatterplot-Location: Langlev-Element: Ni. 
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Figure A.17 Graphite Tube Scatterplot-Location: Langley-Element:S1. 
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Figure A.{8 Graphite Tube Scatterplot-location: Langley-Element: Ti. 
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Figure A.19 Graphite Tube Scatterplot-Location: Myrtle Beach-Element: le. 
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Figure A.20 Graphite Tube Scatterplot-Location: Myrtle Beach-Element:A g. 
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Figure A.21 Graphite Tube Scatterplot-Location: Myrtle Beach-Element:Al. 
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Figure A. Graphite Tube Scatterplot-Location:Myrtle Beach-Element:Cr. 
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Figure A.23) Graphite Tube Scatterplot-Location: Myrtle Beach-Element:Cu. 
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Figure A.24 Graphite lube Scatterplot-Location: Myrtle Beach-Element: Mg. 
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Figure A.25 Graphite Tube Scatterplot-Location: Myrtle Beach-Element:Ni. 
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Figure A.26 Graphite Tube Scatterplot-Location: Myrtle Beach-Element:Si. 
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Figure A207) Graphite Tube Scatterplot-Location: Myrtle Beach-Element:Ti. 
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Figure A.28 Graphite Tube Scatterplot-Location:Pensacola-Element: Fe. 
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Figure A.29 Graphite Tube Scatterplot-Location:Pensacola-Element:\g. 
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igure A.30) Graphite Tube Scatterplot-Location:Pensacola-Element:Al. 
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Figure A.31 Grapiite Tube Scatterplot-Location:Pensacola-Element:Cr. 
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Figure A.32 Graphite Tube Scatterplot-Location:Pensacola-Element:Cu. 
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Figure A.33 Graphite Tube Scatterplot-Location:Pensacola-Element: Mg. 
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Figure A.34 Graphite Tube Scatterplot-Location:Pensacola-Element:N1. 
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Figure A.35 Graphite Tube Scatterplot-Location:Pensacola-Element:Si. 
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Figure A.36 Graphite Tube Scatterplot-Location:Pensacola-Element:Ti. 


67 


APPENDIN B 


ANOVA TABLES FOR CORRELATION DATA 
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APPENDIX D 
REGRESSION MODEL 


The primary regression model used in this paper was the siunple linear regression 
model with least squares estimators. The equation vielded from this model will give the 
best linear unbiased estimates for the data when all assumptions are met. The general 


Straight line regression model is given in Equation D.1. 


Vo = 0 7 ee (eqn 1.1) 
The assumptions used for this model are [Ref. Il: p. 465): 
1) Hee have a population of y values for each x; the population variable 


corresponding lox, is V,. 


2) iD | = a geil each x; Glas 0). 
3 Marly] = 6° for cach x; (homoscedascity’). 
4) the errors of observation, e; = ¥;- a+ bx; are uncorrelated. 
Solutions to the regression model are obtained through least squares estimation 


techniques. The resultant equations that provide the estimutes are: 


% = AR tat). (eqi 12s) 
MO =) 

a ae ee 

a=) - cae (eqn D.3) 
where 

y = (Lym, (eqn D.4) 
and 

fs (yx )in (eqn D.5) 


‘The variance estimates are: 


Wee 


a = ° 
NE Se 
ey san) 
and 
+ 7 
5 g- Vix. 
a aT: 
yas, - 5 
where 


1S 


(eqn 1.6) 


CSOD 


(eqn D.3S) 
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